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5h ; ABSTRACT 

The search for new phenomena at hadron colliders requires a good understanding 
of QCD processes. The analysis of multi-jet signatures in the top quark search at 
the Tevatron is one example, forward jet-tagging and rapidity gap techniques in 
the analysis of weak boson scattering events at the LHC will be another important 
application. These topics are discussed in the context of multi-parton/multi-jet 
QCD processes. Also described are some of the calculational tools, like amplitude 
techniques and automatic code generation for tree level processes. 
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(TV !• Introduction 

o ' 

Ever since the start-up of the SppS at CERN, hadron colliders have represented 
the high energy frontier in particle physics. Up to this day only they have been able to 
produce lU-bosons at a meaningful rate, even though this situation will change in the 
summer of 1996 with the running of LEP at an energy of 161 GeV. Top quark physics 
is another area which will only be accessible to hadron colliders in the forseeable future, 
namely to pp collisions at the Tevatron for the next decade and to pp collisions at the 
LHC afterwards. Even the construction of an e + e~ linear collider or a collider 
will leave the LHC as the machine with the highest usable center of mass energy and, 
hence, at the forefront of discovery physics. 

The discovery of W and Z bosons at the SppS0 with their clean leptonic decays 
may well represent the exception rather than the rule for the observation of new physics 
at hadron colliders. The top quark search,!~S in particular in the ti — > bW + bW~ — > 
t^v + 4 jets decay mode, is more representative of the complex signatures which must 
be searched for in order to discover new particles. The production and decay of squarks 
and gluinos, for example, would lead to multi-particle final states which would give rise 
to multiple jets and missing transverse momentum signatures or to final states with 
several leptons and jets in the case of cascade decays of the supersymmetric particles.! 

In order to distinguish these new physics signals from Standard Model (SM) back- 
grounds, typically produced by the emission of extra, hard quarks or gluons in simple 
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electroweak processes, the precise features of these multi-parton production processes 
must be determined, both for the signals and for the backgrounds. Fortunately, am- 
plitude techniques have been developed over the last 15 yearsBu which are very well 
suited for this task and which allow to perform parton level Monte Carlo studies with 
full tree level matrix elements. These lectures describe some recent uses and some of 
the developments in this field of QCD multi-parton production processes. 

In Section the top quark search at the Tevatron will be discussed as an example, 
showing why the study of multi-parton processes is needed: in order to establish the 
signal and to measure the top-quark mass kinematically, the W + 4 jets background, 
with and without b-quarks in the final state, had to be understood. Motivated by this 
example, Section [3] then provides an introduction into the calculation of multi-jet cross 
sections. I start with a basic discussion of parton level Monte Carlo programs and then 
describe a particular amplitude technique, developed by Hagiwara and myself,Q for the 
efficient evaluation of partonic cross sections. This technique has been encoded in the 
HELAS program packages and is being used by MadGraphli3 to automatically generate 
the FORTRAN code for (almost) arbitrary SM transition probabilities. The basic ideas 
behind these programs will be described. 

An important application of multi-parton processes is the study of weak boson 
scattering at the LHC,Hj which includes Higgs production by weak boson fusion, e.g. 
the process qq — > qqH —> qqW + W~ . In Section || the production and subsequent decay 
H — > W + W~ of a heavy Higgs boson will be used as an example to discuss some 
of the techniques, like forward jet-tagging&ll! or central jet vetoing, EE which have 
been developed to distinguish weak boson scattering signals from QCD background 
processes. 

Of particular interest here are the different color structures of weak boson scat- 
tering (no color exchange in the t-channel) as compared to typical backgrounds where 
color is exchanged between the two incident partons. This different color structure 
leads to very different patterns of radiated gluons.00 In the signal, gluons are emitted 
at rather low transverse momenta and mainly in the forward and backward direc- 
tions, leaving a central region with little hadronic activity apart from the Higgs decay 
products. Typical background processes emit gluons at considerably higher transverse 
momenta and preferentially in the central region. These differences suggest techniques 
like a minijet vetdia or rapidity gapstll to enhance the signal versus the backgrounds. 
They will be studied in Section [| It is helpful, however, to first consider gluon radiation 
patterns and the consequences of t-channel color singlet exchange in a simpler environ- 
ment, namely dijet events at theTevatron. These topics and the recent observation of 
rapidity gaps at the Tevatrontacll will be discussed in Section f|. 

2. Jets in the Top Quark Search 

The discovery of the top quark§@ at the Fermilab Tevatron provides a beautiful 
example for the use of hadronic jets as a tool for discovering new particles. Let us have 
a brief look at the top quark search at the Tevatron from this particular viewpoint: 
the use of multi-parton cross section calculations as a necessary ingredient in particle 
searches. 
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In pp collisions at the Tevatron, at a center of mass energy of 1.8 TeV, the top 
quark is produced via quark anti-quark annihilation, qq — > ti, and, less importantly, 
via gg — > ti. Production cross sections have been calculated at next-to-leading orderB 
and are expected to be around 4 pb for a top mass of 180 GeV.il The large top decay 
width which is expected in the SM, 

r(t -> W + b) « 1.7 GeV , (1) 

implies that the t and t decay well before hadronization, and the same is true for the 
subsequent decay of the W bosons. Thus, a parton level simulation for the complete 
decay chain, including final parton correlations, is a reliable means of predicting detailed 
properties of the signal. In order to distinguish the top quark signal, ti — > bW + bW~ , 
from multi-jet backgrounds, the leptonic decay W — > tv {t = e, /i) of at least one of 
the two final state Ws is required. Since the leptonic decay of both Ws has a branching 
ratio of ~ 4% only, the prime top search channel is the decay chain 

ti -> bW + bW- -> tvqqbb, (2) 

which, within the SM, has an expected branching ratio of ~ 30%. After hadronization 
each of the final state quarks in (H) may emerge as a hadronic jet, provided it carries 
enough energy. Thus the ti signal is expected in W + 3 jet and W + 4 jet events^} 

Events with leptonic W decays and several jets can also arise from QCD correc- 
tions to the basic Drell-Yan process qq — > W — > I v. The process ug — > dggW + , for 
example, will give rise to W + 3 jet events and its cross section and the cross sections 
for all other subprocesses with a W and three partons in the final state need to be 
calculated in order to assess the QCD background of W + 3 jet events, at tree level. 
W + n jet cross sections have been calculated for n = 3 jetso and n = 4 jets! and 
some of the methods used in these calculations will be discussed in Section [3|. As in the 
experiment, the calculated W + n jet cross sections depend critically on the minimal 
transverse energy of a jet. CDF, for example, requires a cluster of hadrons to carry 
Et > 15 GeV to be identified as a jet,@ and this observed Et must then be translated 
into the corresponding parton transverse momentum in order to get a prediction for 
the W + n jet cross sections. 

At this level the QCD backgrounds are still too large to give a viable top quark 
signal. The situation is improved substantially by using the fact that two of the four 
final state partons in the signal are 6-quarks, while only a small fraction of the W + 
n parton background events have 6-quarks in the final state. These fractions are readily 
calculated by using W + n jet Monte Carlo programs. There are several experimental 
techniques to identify 6-quark jets, all based on the weak decays of the produced 6's. 
One method is to use the finite b lifetime of about r = 1.5 ps which leads to 6-decay 
vertices which are displaced by 7cr = few mm from the primary interaction vertex. 
These displaced vertices can be resolved by precision tracking, with the aid of their 

*Gluon bremsstrahlung may increase the number of jets further and thus all W+ > 3 jet events are 
potential ti candidates. 
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Fig. 1. Number of W + n jet events in the CDF top quark search as a function of jet 
multiplicity. Number of observed events are given without 6-tagging (open circles) and with 
an SVX tag (triangles). The expected background, mainly from QCD W + n jet events, is 
given by the cross-hatched bars. From Ref. 3. 



Silicon VerteX detector in the case of CDF, and the method is, therefore, called SVX 
tag. In a second method, b decays are identified by the soft leptons which arise in the 
weak decay chain b — > W*c, c — > W*s, where either one of the virtual Ws may decay 
leptonically.™ 

The combined results of using jet multiplicities and SVX 6-tagging to isolate the 
top quark signal are shown in Fig. |I} A clear excess of 6-tagged 3 and 4 jet events is 
observed above the expected background. The excess events would become insignificant 
if all jet multiplicities were combined or if no 6-tag were used (see open circles). Thus 
jet counting and the identification of 6-quark jets have been critical for the discovery 
of the top quark. 

Beyond counting the number of jets above a certain transverse energy, the more 
detailed kinematic distributions, their summed scalar and multi-jet invariant 

masses, have also been critical in the top quark search and thus needed to be predicted 
for the W + n jet backgrounds. The top quark mass determination, for example, relies 
on a good understanding of these distributions. Ideally, in a tt — > (£ + vb)(qqb) event, for 
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Fig. 2. Reconstructed mass distribution for W+ > 4 jet events with a 6-tag. The solid 
histogram represents the CDF data. Also included are the expected background (dotted 
histogram) and the expected signal+background for m t = 175 GeV. The insert shows the 
likelihood fit to determine the top quark mass, which yielded mt = 176 db 8 db 10 GeV. From 
Ref. 3. 



example, the two subsystem invariant masses should be equal to the top quark mass, 

m t ~ m(l + vb) « m(qqb) . (3) 

Including measurement errors, wrong assignment of observed jets to the two clusters, 
etc. one needs to perform a constrained fit to extract m t . The CDF result of this fitfl 
is shown in Fig. ^. Combining this with the corresponding DO measurements gives the 
valued 

m t = 180 ± 13 GeV . (4) 

In addition, Fig. || demonstrates that the observed 6-tagged W + 4 jet events (solid 
histogram) are considerably harder than the QCD background (dotted histogram). 
On the other hand the data agree very well with the top quark hypothesis (dashed 
histogram). The kinematic distributions of the jets contain much valuable information 
and in order to extract it we must be able to reliably calculate them for processes with 
many partons in the final state. 
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3. Monte Carlos and Amplitude Techniques 

In the previous Section we have seen that the isolation of new physics signals 
like top quark pair production requires precise predictions for the properties of signal 
and background processes. One needs to know how the corresponding differential cross 
sections change as a function of the kinematic variables of the observed leptons and 
jets. Starting with a generic process 

pp -> B 1 B 2 . . . B n X , (5) 

where the Bi stand for the final state leptons and jets, we will discuss in the follow- 
ing how to determine the corresponding cross sections in a parton level Monte Carlo 
program, at tree level. We begin by considering the general structure of a Monte Carlo 
program in Section 3.1. A formalism for the fast numerical evaluation of polarization 
amplitudes is described in Section 3.2, and Section 3.3 discusses automatic code genera- 
tion for the calculation of squared matrix elements, as implemented in the MadgraphtJ 
program. 

3.1. Monte Carlo simulations at the parton level 

At the parton level, the process pp — > B\B 2 . . . B n X originates from a number of 
different subprocesses, 

a x a 2 -> hb 2 • • • b n X , (6) 

and the cross sections for all these subprocesses must be added to obtain the probability 
for observing the final state B\B 2 ■ ■ ■ B n . 

Consider the B\B 2 = W + j signature as an example. The contributing parton 
level subprocesses a\a 2 — > bib 2 are 



du, sc, ud, cs — > W + g , (7) 

ug, gu -> W + d , (8) 

eg, gc -> W + s , (9) 

dg, gd -> W + u , (10) 

sg, gs - W + c, (11) 



where Cabibbo mixing has been neglected for simplicity^]. The various subprocesses 
are closely related, of course: by choosing either the proton or the anti-proton as the 
source of the initial state quark or anti-quark, by interchanging first and second gener- 
ation quark flavors, or by crossing the gluon between the initial and final states. Using 
these relations, the required analytical calculations can be reduced significantly. At the 
numerical level, however, only the flavor independence of QCD corrections allows to 
reduce the number of independent cross section evaluations. 

^Because of the small value of the Cabibbo angle, sin 2 ^ = 0.049, the resulting error is negligible 
compared to the neglect of higher order QCD corrections, for example. It is further mitigated by the 
fact that after summing over final state flavors and neglecting quark masses, only the annihilation 
processes (R) depend on the value of the Cabibbo angle. 
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In order to calculate the cross section for the BiB 2 . . . B n final state, the squared 
matrix elements, summed over final state colors and polarizations and averaged over 
initial ones, 

El^r = ^ color L a) E E IMI 2 (12) 

TT \ 1 2 J co l ors polarizations 

must be known for all contributing subprocesses. Their determination is the subject of 
the next two subsections. The full cross section is then given by 

a = I d,Xidx 2 fax/pi^l) fa 2 /p(x 2 ) 

subprocesses 

d$ n (xip + x 2 p; pi ■ ■ .p n )e(cuts) ^2\M\ 2 (aia 2 -> hb 2 • • • b n ) . (13) 

Here f ai / p (xi) is the probability to find parton a\ inside the proton, carrying a fraction 
x\ of the proton momentum, i.e. the ai momentum is X\p. Similarly, f a2 /p( x 2) is the 
parton a 2 distribution function inside the anti-proton. In the second line of (|13|) l/2s 
is the flux factor for the partonic cross section, and 

rf$ n (P; p x . . .p n ) = ft ((^%) (2^) 4 ~ Eft) (14) 

is the Lorentz invariant phase space element. Finally, G(cuts) is the acceptance func- 
tion, which summarizes the kinematical cuts on all the final state particles, i.e. = 1 
if all the partons a, satisfy all the acceptance cuts and = otherwise. 

The calculation of the cross section in eq. (|13|) involves a 3n — 2 dimensional 
integral which is best performed by Monte Carlo. In order to do this one first maps the 
integration region onto a 3n — 2 dimensional hypercube of unit length, i.e. one rewrites 
the integration measure as 

-y 3n-2 

—dxidx 2 d& n = J Yl ' (15) 
2 s i=i 

where J is the Jacobian factor of the transformation of integration variables. The Monte 
Carlo integration then samples the integrand randomly at some large number, N, of 
points in the hypercube. In other words, take iV sets {rj} of 3n — 2 random numbers 
each. A good approximation to the cross section is then provided by 

cr^i^J £ /(*0/(* 2 )E|.M| 2 e(cuts) , (16) 

{n } subprocesses 

which approaches the correct cross section value in the limit — > 00. 
Eq. ( |16|) contains all the elements of a Monte Carlo program: 

• An integration routine (for example the programs VEGAS@ or BASESEa) pro- 
vides the set {rj} of random numbers and the weight factor wgt = U l/N". Each 
set {rj} is called an event. 
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• The phase space generator generates the kinematic variables, i.e. the parton mo- 
menta Pi, the initial parton momentum fractions x\ and x 2 , and the Jacobian 
factor J from this set of random numbers. 

As a trivial example consider the two body phase space in the decay process a(P) — > 
b\(pi)b 2 (p 2 ) , where b\ and b 2 represent two massless particles. The Lorentz invariant 
phase space element is best written in polar coordinates in the P-rest frame: 

d$ 2 (P; Pl ,p 2 ) = -L d n = ^12 (it) 

where the components of the 4-momentum pi are given by the two random numbers 
r\ and r 2 via 

VP* 

Pi = — - — (1, sin 9 cos 0, sin 9 sin 0, cos 9) , (18) 
cos 9 = 2n-l , (f) = 2vrr 2 . (19) 

Since the phase space element rf$ 2 is Lorentz invariant, the momenta may now be 
boosted to an arbitrary frame and the momentum p 2 is then calculated from momentum 
conservation, 

p 2 = P- Pl . (20) 

• Given the 4- momenta of the final state particles {p{\ one can easily check if 
the event passes the acceptance cuts, i.e. 0(cuts) = 1. If the event is accepted 
then the squared matrix elements, Sl-M| 2 > are calculated for all contributing 
subprocesses. 

• The weight factor (value of the integrand) 

w = J £ f(x 1 )f(x 2 )Y / \M\ 2 e(cuts) (21) 

subprocesses 

is returned to the integration routine, which then adds these weights for the 
various regions of the hypercube, thus obtaining an estimate of the cross section, 
a, according to eq. (p~6|) . 

• At the same time one can generate arbitrary distributions da/dz by filling a 
number of histograms. For each event one uses the 4-momenta pi to determine 
the value of z for each histogram in question. Given the bin width Az one adds 
w -wgt/Az to the bin corresponding to z and thus obtains a correctly normalized 
histogram of da/dz at the end. 

In the above procedure the choice of mapping from the unit hypercube to the 
Lorentz invariant phase space element is obviously ambiguous. The art of writing a 
good phase space generator is to choose a mapping which results in a large fraction of 
events having a sizable weight, close to the maximum weight w max encountered in the 
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full sample. Events with small weights contribute very little to a in eq. (|16D and thus 
do not reduce the statistical error, but they require as much time to be computed as 
high weight events. The task in writing the phase space generator is thus to anticipate 
the dominant phase space regions, those regions in which poles in the matrix elements 
lead to large cross sections, and to find a mapping which smooths out these maxima in 
the product JJ2\M\ 2 . It is a good practice to histogram the distribution of the weight 
factors, w, in order to judge whether improvements to the phase space generator should 
be made. 

3.2. Amplitude techniques 

It is the squared amplitudes, D|-M| 2 , which contain the full information on the 
dynamics of the underlying physical processes and which are at the core of any cross 
section calculation. For each subprocess a%(ki) a 2 (&2) —* &i(pi) ■ ■ • b n (p n ) and phase 
space point (ki, k 2 ; Pi, ■ ■ ■ , p n ) the squared amplitude has a well defined numerical value. 
The question then is: What is the most efficient way to "evaluate the Feynman graphs" 
and to thus obtain this value? 

Traditionally, the method of choice has been to use "trace techniques" to express 
X)|.M| 2 in terms of scalar products p, L ■ pj. For relatively simple processes, like 2^2 
scattering, one thus obtains very compact and easy to interpret analytical expressions, 
in particular for unpolarized scattering, and no residual polarization and color sums 
need to be performed in evaluating Eq. ([12]) numerically. Unfortunately, the expressions 
become very large as the number, /, of individual Feynman graphs increases. For 

/ 

M = Y t M i (22) 
1=1 

a total of /(/ + l)/2 cross terms Re(A4*J\Aj) need to be calculated and hence the 
complexity of the problem grows at least with the square of the number of Feynman 
graphs. For 2 — > 4 or 2 — > 5 processes the task typically becomes forbidding. 

For such complex processes the resulting analytical expressions cannot easily be 
interpreted any more, and in practice one merely needs them for subsequent numer- 
ical evaluation. It then becomes advantageous to numerically evaluate the individual 
Feynman amplitudes Mi as / complex numbers which can be summed and squared 
trivially. The complexity of this problem only grows linearly with /. 

A number of approaches have been developed for the direct evaluation of polar- 
ization amplitudes. The CALKUL method! and improvements on itS may be the most 
widely used and it yields rather compact expressions for helicity amplitudes involving 
massless fermions. Similar (high) numerical speeds are achieved via the HZ method^ 
which was developed with massive fermions in mind, but which simplifies significantly 
as well when massless fermions are considered. As compared to CALKUL results, he- 
licity amplitudes in the HZ method are typically not as compact. However, the basic 
structure of the underlying Feynman graphs is preserved and this fact has proven in- 
valuable in the automation of code generation.^ At the same time it becomes very easy 
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to add new features to existing calculations, like implementing finite width effectsc3 or 
anomalous gauge boson couplings.0 

The HZ method is implemented in a FORTRAN77 package called HELAS.tH In 
the following the principles behind the method will be presented. For full details the 
reader is referred to the original articles.iil 

In order to exploit the simplifications which occur for massless fermions we work 
in the chiral representation, with 75 given by 

*- (";:)■ «> 

External spinors are split into their chiral components, 

*-(£)■ (24) 

where ip- and ip+ are two-component Weyl spinors of negative and positive helicity, 
respectively. As an example consider the wave function of an incoming fermion of 4- 
momentum p^ = (p° ,p x ,P y ,Pz) = OA p) and helicity a/2. It can be represented by a 
M-spinor with chiral components 



u(p, a)± = \fp° ± a\p\ Xa(p) • (25) 
Here the Weyl spinors of fixed helicity a/2 = ±1/2 are given explicitly by 

'2 p p +p z ) \ Px + Wy J 



- w ■ ^rra ( " « - ) ■ 

Given this explicit representation, the four components of the w-spinor can be calculated 
as a set of four complex numbers once the external 4-momentum and the helicity are 
fixed. In the HELAS package this is achieved by a simple subroutine call, 

call IXXXXX(P,FMASS,NH,+1,PSI) 
where P=P(0:3) is the external fermion 4-momentum, FMASS is the fermion mass, 
NH= a = ±1 its helicity, the entry +1 indicates that an external fermion (as opposed 
to anti-fermion) wave function needs to be calculated and the output PSI = PSI(1 :6) 
contains the resulting spinor ip = u(p, a) and the 4-momentum p stored as 

PSI = (V>i, 1P2, fa, ^4, P° + iPz, Px + ip y ) ■ (28) 

In the same way the wave function v(p, a) of an anti-fermion in the initial state is 
calculated via 

call OXXXXX (P, FMASS, NH=cr,—l ,PSI) 
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where the entry —1 indicates that the wave function of an anti-fermion is needed and 
the name OXXXXX refers to the direction of the arrow on the fermion line in the Feynman 
graph: an outgoing fermion line is considered. 

The other wave functions which are needed in the calculation of Feynman ampli- 
tudes are the polarization vectors of external vector bosons. Given the 4-momentum 
k^ 1 = (E, k x , k y , k z ) = (E, k) of a vector boson of mass m, three independent polar- 
ization vectors which satisfy k ■ e=0 are given by (k? = \Jk"l + ky) 

e"(A;, 0) = -^( k V^, k x , k y , k z ) , (29) 
m\k\ 

£^{k, 1) = jj^-j-^ (0, k x k z , kyk z , k T ) , (30) 

= -^(0, -ky, k x , 0) . (31) 

The polarization vectors for fixed helicities A = ±1 are obtained from this Cartesian 
basis via 

e»(k, \ = ±l) = -L ( T£ »(k, 1) - &(k, 2)) (32) 

and they are obtained in HELAS via subroutine calls 

call VXXXXX(K,VMASS,NHEL(= A),-1,EPS) 
where —1 indicates that the wave function for a vector boson in the initial state is to 
be calculated (+1 would be appropriate for calculating the wave function e^ 1 * of a final 
state vector boson) and as for fermions EPS=EPS(1:6) combines the four components 
of the complex polarization vector with momentum information on the vector boson. 

Having determined the external particle wave functions we need to do the 7- 
matrix algebra corresponding to the fermion lines inside Feynman graphs. For any 
4-vector a M the contraction with 7^, in the chiral representation, can be written as 

* = { ,1 (33) 



Here the 2x2 submatrices {<f)± are given by 

-r( n — t„ \ \ 

(34) 



a°Ta z ^(a x -ia y ) 



\ T(a x + ia y ) a ± a z 

Now consider external vector bosons attached to a fermion line as depicted in 
Fig. ^. The external fermion spinor ip = ip(p), the first vector boson emission and the 
fermion propagator of momentum p — k\ give rise to a new four-component spinor, 

\k\, V >= -j 7T h (-Wl n 75 - 19 r + r) ry5 ) i> 5 ( 35 ) 

f — f x — rrif V 2 2/ 
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Fig. 3. Piece of a Feynman graph corresponding to vector boson emission off a fermion line. 



where the left- and right-handed couplings gt and g^ depend on the nature of the 
vector boson considered. In our explicit representation of Dirac-matrices and spinors 
the resulting ket is obtained by simple 4x4 matrix multiplication which is performed 
numerically in HELAS via a subroutine call 

call FVIXXX (PSI , EPS1 , G , FMASS , , PSIV) 
Here PSI and EPS1 contain the full information on the external spinor and the vector 
boson of momentum k%, G = (—<?£, —Qr) stands for the coupling constants, the 
indicates that the fermion decay width has been set to zero and the output PSIV = 
PSIV(1 : 6) again stores the four-component spinor \ki, p > and the 4-momentum p—k\ 
on the outgoing spinor line, encoded in terms of two complex numbers. 

The emission of additional vector bosons now becomes trivial: the output PSIV 
of the last call of FVIXXX is used as the input PSI of the next FVIXXX call, with vector 
boson 1 replaced by vector boson 2, etc. This leaves us with the last vector boson 
vertex, 

V>(P2)^ n (~Wl — - igR ~^ 75 ) I* 1 ' • • • ' fc"- 1 ' p > ' ( 36 ) 

i.e. the task of calculating the final amplitude, corresponding to the Feynman graph 
of Fig. |3|, from the complex 4-spinor \k±, . . . , fc n -i> p >, the complex conjugate spinor 
ip{p2) and a 4x4 matrix sandwiched in between. This, of course, is a trivial task and 
accomplished in HELAS by the call of another subroutine, IOVXXX. Alternatively, if 
the last vector boson is virtual, one may want to calculate the 4-vector 

- my \ m v J V 2 2 / 

(37) 

which can then be used in place of the polarization vector for an external vector boson 
in e.g. the calculation of a second fermion line. This is done by the subroutine JIOXXX. 

In an analogous way non-abelian couplings of vector bosons are handled. From the 
polarization vectors and momenta of two external gauge bosons (or the corresponding 
currents in the case of virtual gauge bosons) one may calculate a new current (call 
JVVXXX) . Alternatively, one can calculate an amplitude from the "wave functions" and 
momenta of three gauge bosons, via the call of yet another subroutine, VVVXXX. 

In this fashion numerical expressions for arbitrary Feynman diagrams can be 
built in a block-like fashion. This method actually is quite economic, since many of the 
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Fig. 4. Feynman graphs for the process e + e — > uugg. 

building-blocks will be common to several Feynman graphs and only need to be calcu- 
lated once. As an example consider the Feynman graphs for the process e + e _ — > uugg 
depicted in Fig. |]. The incoming e + e~ current is common to all 8 graphs. Attaching 
this current to the external w-spinor of the u leads to a 4-spinor which is common to 
graphs 1, 2, and 3. The virtual gluon "current", which describes the 3-gluon vertex, 
appears in graphs 2 and 5, and so on. Exploiting the repetitive elements in the various 
Feynman graphs substantially increases the efficiency of the numerical calculation, in 
particular for processes with a very large number of graphs. 

Finally one should mention another convenient feature of the HELAS program. 
When crossing initial and final state particles the fermion arrows in the Feynman graphs 
do not change, only the overall signs of the 4-momenta of the crossed particles need to 
be flipped. This is achieved in HELAS by flipping the sign factors in the determination 
of the external wave functions: In the IXXXXX, OXXXXX and VXXXXX calls discussed 
above the sign of the second last argument needs to be reversed. No further changes 
are required when calculating crossing related processes. 

3.3. Automatic code generation: the MadGraph example 

The discussion of the previous subsection should have made it clear that writing 
the code, which evaluates the squared amplitudes for the various subprocesses, is actu- 
ally quite a mechanical and boring exercise, exactly the type of task which should be 
relegated to a computer. Indeed, a number of programs have been developed in recent 
years for automatic code generation for scattering processes.0H The example to be 



13 



discussed below is MadGraph.tO MadGraph is a standard Fortran program^ which gen- 
erates Fortran code calling the HELAS subroutines. At present MadGraph is capable of 
generating the code for the squared matrix elements for SM processes with up to seven 
external particles. It already has been used successfully in state of the art calculations 
like Z + 4 jets production in pp collisions. 




Fig. 5. Allowed topologies for Feynman graphs with four external particles. The circles 
indicate the 25 locations where a fifth external particle could be attached. 



Let us consider the various steps which MadGraph takes to generate a program 
for the calculation of I]|-M| 2 - The first step in calculating the scattering amplitude M. 
is to generate all possible topologies for tree level Feynman graphs, keeping in mind 
that the SM, as a renormalizable field theory, allows 3 and 4 particle vertices only. 
As an example, the four possible topologies for Feynman graphs with four external 
particles are shown in Fig. [5[ A fifth external particle could be attached to each of 
these graphs in the locations indicated by the open circles, leading to 25 distinct tree 
level topologies for Feynman graphs with 5 external particles. For larger numbers of 
particles the number of allowed topologies grows very rapidly, as shown in Table |l|. This 
rapid growth leads to increased memory demands of MadGraph and puts a practical 
limit of 7 external particles on present day workstations. 



Table 1. Number of distinct topologie of tree level Feynman graphs in a renormalizable 
field theory as a function of the number of external particles. 



# external particles 


3 


4 


5 


6 


7 


8 


# topologies 


1 


4 


25 


220 


2485 


34300 



The tree level amplitude is uniquely fixed once the external particles and the 
interaction lagrangian are specified^]. In general only a subset of the allowed topolo- 
gies will actually contribute to a given subprocess. Only the third topology in Fig. |5| 
contributes to e + e~ — > for example, but with two distinct Feynman graphs, cor- 

responding to photon and Z exchange. MadGraph contains a table of all allowed SM 

* The MadGraphEl source code and the HELAsl routines called by the Madgraph output can 
be obtained via anonymous ftp from phenom.physics.wisc.edu where it is located in the directory 
pub /madgraph/ . 

§ For SM subprocesses involving four or more external quarks both gluon and electroweak boson 
exchange is possible between the quark lines. In these cases, MadGraph allows the user to choose the 
desired order in the strong coupling constant a s . 
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vertices. For a given topology it first assigns the external particles to the external lines 
in all permutations. Madgraph then works inwards and checks at each vertex with one 
unknown line whether there is an allowed assignment of the not yet labeled internal 
propagator, according to the vertex table. Storing the allowed assignments, the process 
is continued until either no allowed particle assignment for the last propagator can 
be found at some vertex or all propagators have been labeled. In the latter case the 
complete assignment is stored as a Feynman graph which needs to be calculated. 

The HELAS code generation follows a similar scheme. First the external particle 
wave functions are calculated (i.e. the appropriate calls of IXXXXX, OXXXXX or VXXXXX 
are printed in the output program). The wave functions are, of course, common to 
all the Feynman graphs. For a vertex with one unknown line, one of the function 
calls described in the previous subsection then determines the (virtual particle) wave 
function corresponding to this line, be it a current for a gauge boson or a 4-spinor for a 
fermion. This procedure ends at a last vertex where all wave functions of the incident 
particles have been calculated. One final subroutine call then determines the complex 
value, Ai, of the amplitude for the given Feynman graph. 

The Ai only contain the propagator factors, coupling constants and the Lorentz 
structure of the vertices. For QCD processes the color factors are not yet included. 
For a given Feynman graph the color factor is a product of factors from each vertex^. 
These are, for example, 





vertex : 


\a 

ij 2 ' 




(39) 


9 a 9 b 9 C 


vertex : 


• fabc ) 




(40) 


g a g b g c gd 


vertex : 


jabe jcde 


(3 vertices) , 


(41) 


Qi Qj 7 


vertex : 






(42) 



In the squared amplitude all colors are summed/averaged over and MadGraph uses the 
familiar commutator and completeness relations 

% f abc = 2tiT a [T\T c ] , (43) 

T$I% = \5 a 5 k] - ^5 l3 5 kl , (44) 

to calculate the color summed/averaged products, C^-, of the color factors for any pair 
of Feynman graphs, Thus, the squared amplitude takes the form 

Ei-^i 2 = E (45) 

^Each 4-gluon vertex leads to three distinct Feynman graphs in MadGraph, corresponding to the 
three different color structures in the Feynman rule 

- f abe f cde (g w 9u, - g^g vp ) - f ace f dbe {g^g vp - g^g« P ) - f ade f bce (g^g P * - g w g va ) ■ (38) 
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A final trick is then used to simplify the color structures: the hermitian matrix is 
diagonalized, yielding n max nonzero eigenvalues A n and orthogonal eigenvectors z n j, 

f 

^2 CijZnj — Z* n i ■ (46) 

i=i 

This allows to write the squared amplitude in terms of a much smaller n max x / 
dimensional matrix z and a set of n max eigenvalues A n , 

Tlmax f 

£|.M| 2 = £ A n |£ z ni Atf . (47) 

n=l i=l 

In MadGraph the eigenvalues A n and the z n i are also used to absorb additional 
factors which appear in the cross section formulas for specific processes: 

• A relative factor —1 between subamplitudes which are obtained by interchanging 
external fermions. 

• A factor 1/JV! for each set of N identical particles in the final state. 

• An overall factor initial colors)" 1 from taking the average over initial state 
colors. 

These last two factors must be kept in mind when using the simple crossing relations of 
the HELAS functions to obtain squared amplitudes for crossing related processes: Mad- 
Graph generates the code for calculating the polarization and color averaged amplitude 
squared for a specific input process. 

Using MadGraph actually is very simple, one merely needs to type in the desired 
process in a self-explanatory shorthand fashion. It is prudent, however, to keep the 
general MadGraph and HELAS properties in mind, so as to be able to change the 
generated programs to ones needs. Desired changes may involve the addition of new 
physics contributions to the amplitudes, a more sophisticated treatment of finite width 
effects, or a more economic programming of crossing related processes. One should 
always keep in mind that automated systems are merely intended to free us from the 
labor of mindless coding. They are tools which make studying the dynamics of multi- 
parton processes more enjoyable and which allow us to concentrate on all the physics 
aspects of signal and background processes. 

4. Manifestations of Color: Rapidity Gaps at the Tevatron 

One important aspect distinguishing signal and background processes is the color 
flow in the contributing Feynman graphs. Many signal processes involve electroweak 
interactions and are thus due to the exchange of color singlet quanta while the QCD 
backgrounds are mediated by color octet gluon exchange. The question arises whether 
this different color flow has observable consequences which can then be used for back- 
ground reduction. There is no general answer to this question. However, there is an 
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important class of processes, mediated by t-channel color singlet exchange, where the 
answer is a clear "yes". One of the most important fields of study at the LHC, weak 
boson scattering, belongs to this class and will be addressed in more detail in the next 
Section. Let us start here with a much simpler example, namely dijet production at 
the Tevatron, which has emerged as an unexpectedly rich QCD laboratory in recent 
years. The discussion will closely follow Ref. 35. 



h>Pi 



:y,z 



Q Q 

(a) 



i 2 ,p 2 



i 4 ,P 4 



h>Pi 



l y P 3 



8 



Q Q 
(b) 



i 2 ,p 2 



i 4 ,P 4 



Fig. 6. Feynman graphs contributing to elastic quark-quark scattering via a) color singlet 
photon or Z exchange and b) color octet gluon exchange. 



4-1. Quark-quark scattering at the tree level 

A number of 2 — > 2 subprocesses contribute to dijet production: gluon-gluon 
scattering, quark-gluon scattering, qq annihilation and pair production, and quark- 
quark or quark- ant iquark elastic scattering. For simplicity let us concentrate on the 
Rutherford scattering process 

Ql(jPl,il)Qs{P3,i3) -» Q2(P2,i2)Q4(P4,h) , (48) 

where the i n and p n denote the colors and momenta of the quarks. The two types of 
SM contributions, namely color singlet 7/Z exchange and color octet gluon exchange 
in the t-channel are depicted in Fig. |6|. Neglecting the Z contribution in the following, 
the matrix elements for the two cases can be written as 

KA R t\ p p 4 - A A pp ^b^bi) ^{p^^u^) 

(P1-P2) 2 

M g = T? 2il T? 4i3 g 2 A= <^5 i2i3 6 i4il - ^<Ai.) 9 2 A . (50) 

The Kronecker deltas in M. 1 describe the fact that the color of q\ is directly transferred 
to q2 while in the QCD case the exchanged gluon transmits the color of qi to Q 4 
(to leading order in 1/N). Does this difference in color structure lead to observable 
consequences? 

A qualitative understanding of these differences can be obtained by considering 
the acceleration of the color charges in forward scattering, at small \t\ = —{pi — P2) 2 — 
s/2 (1 — cos 6) « (pi + p 3 ) 2 = s. In the case of t-channel photon exchange, quark 
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q will be deflected by a small angle 9 with respect to the original q\ direction and 
so will the color charge carried by it. The acceleration of the color charge by a small 
angle results in gluon bremsstrahlung in the forward direction only and, similarly, 
the small deflection of the color charge carried by Q results in gluon radiation in the 
backward direction, with very little radiation between the two scattered quarks. For 
color octet gluon exchange, on the other hand, the color charges are exchanged between 
the two quarks q and Q, i.e. they are deflected by the large angle n — 9, and this strong 
acceleration of color charges results in strong gluon bremsstrahlung at large angles, 
between the directions of the final state quarks. 

g 




<?i » ' » 5 — * — i 2 



G, — - — 2 — - — Q A 




A m A (2) 



% — *■ 3 * <Z 2 

Q, — *" &» <2 4 



A (5) 




A (3) A (4) 
Fig. 7. Feynman graphs for the process qQ — > gQ<? at tree level. 

This qualitative picture is confirmed by calculating the gluon emission corrections 
to qQ elastic scattering, i.e. by considering the process 

q(pi,h) Q(Ps, is) — >■ q(j>2, 12) Q(P4, u) g(k, a) , (51) 

where k and a denote the gluon momentum and color. At any order of perturbation 
theory, the amplitude for this process can be written in terms of four orthogonal color 
tensors, which we denote by Oi, 2 , Si 2 , and S 3 4, 

M = O x M x + 2 M 2 + S 12 M 12 + S 3i M M . (52) 

In an SU (N) gauge theory they are given explicitly by 

S12 = T°- ih S i2h (53) 
^34 = T» 2ii 8 iAi3 (54) 



-2 
N 

O2 = Tl n 5 l2l3 -T^J un . (56) 



Oi = —(Su + S^+T^J^+T^J^ (55) 
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Si 2 marks the process where the quark q keeps its color. Similarly, S 34 multiplies the 
amplitude for t-channel color singlet exchange as viewed from quark Q. Within QCD, 
0\ and O2 correspond to t-channel color octet exchange as viewed from either of the 
two scattering quarks. 

Let us first apply this color decomposition to the tree level QED and QCD ampli- 
tudes. The five Feynman graphs for the QCD process are shown in Fig. [7|. Lumping the 
momentum and helicity dependence of the individual Feynman diagrams into reduced 
amplitudes ■ ■ ■ one obtains for the QCD amplitudes at tree level 

m qcd = V ( A m + A m^^_ A m (57 ) 

M QCD = V ( A (3) +A (4))^Z|! A (34) (5g) 
,3 



M QCD = _±_ ^(1) + A (2) + A (3) + A Q)j (59) 
M QCD = ^ ^(1) _ A ( 2 ) + A (3) _ A (4) + 2A (5)) = Zl. A M (6Q) 



The non-abelian three-gluon-vertex only contributes to the color octet exchange am- 
plitude M 2 . Both M 2 and Mi vanish for t-channel photon exchange, while the color 
singlet exchange amplitudes are given by 

Mi? D = -ge q e Q (A® + A«) = -ge q e Q ^ (61) 
M 3 Q 4 ED = - Wq (A« + A( 2 )) = -^e g e Q A^ , (62) 

where e g and eQ denote the electric charges of the two quarks. 

The resulting color and polarization summed squared amplitudes, 

N 2 — 1 / N 2 — A \ 

Y t \M?=? L Y 1 N £ |M 12 | 2 +|M 3 4| 2 + 2^^|Mi| 2 + 2|M 2 | 2 , (63) 

" polarizations V / 

are shown in Fig. ||(a) for both the QCD and the QED case. Shown is the dependence of 
Yj I -M 2 1 1 8 on the rapidity of the emitted gluon when all other phase space parameters 
are kept fixed, namely the two quarks are held at px = 30 GeV and pseudorapidities 
7] q = 3 and t]q = — 3 and the gluon transverse momentum is chosen to be pt 9 = 
2 GeV. The details of this choice are irrelevant: forward scattering of the two quarks 
is a sufficient condition to obtain the qualitative radiation pattern of Fig. j|. In the 
QCD case the color octet contributions, via the non-abelian amplitude A^ n£i ', lead to 
enhanced gluon emission in the angular region between the two jets. For t-channel 
photon exchange this region is essentially free of gluons due to color coherence between 
initial and final state gluon radiation ^'0: gluon emission into the central region is 
exponentially suppressed as the rapidity distance from the quarks increases. 

Going beyond the parton level scattering process and considering the actual pp 
scattering event, hadronization of the emitted gluons will result in (soft) hadrons which 
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Fig. 8. Rapidity distribution of emitted gluons in qQ — > qQg scattering for fixed final state 
parton transverse momenta of pTq = 30 GeV and pTg = 2 GeV. The quark rapidities are 
fixed at rj q = ±3 (indicated by the arrows). In part a) results are shown for the sum over all 
color structures for single gluon and for t-channel photon exchange. The M34 terms alone, in 
part b), demonstrate the difference between the QED and the QCD color singlet exchange 
terms. Quark charges are taken as e q = eQ = +2/3e. 



will trace the angular distribution of their parent partons. In t-channel color singlet 
exchange, gluon radiation is severely suppressed in the angular region between the 
quarks, and, thus, a similar pattern of soft hadrons is expected. This leaves a rapid- 
ity gap, i.e. a rapidity region between the two quark jets which is essentially void of 
produced hadrons. 

4-2. Rapidity gaps at the Tevatron 

Rapidity gaps in hard dijet events have indeed been observed at the Tevatron.&0 
The DO Collaboration, for example, has studied events with two hard, forward and 
backward jets with 

E T > 30 GeV , \ Vj \ > 2 , Vl ■ V2 < , (64) 

and has then searched for signs of hadronic activity in the pseudorapidity range be- 
tween the two jets. Inside the region of width Arj c between the tangents to the two jet 
definition cones (of radius R = 0.7 in the lego-plot) DO has then studied the hadronic 
multiplicity in terms of the number of active towers in the electromagnetic calorimeter, 
above a minimum transverse energy threshold of 200 MeV. The measured multiplicity 
distribution is shown in Fig. || One observes a clear excess of low multiplicity events, 
above the expectation from extrapolating a double negative binomial distribution which 
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fits very well the high multiplicity region. This excess corresponds to a fraction 

f gap = 1.07 ± 0.10(stat)+^(syst) % (65) 

of rapidity gap events in dijet production. In addition, this fraction is essentially in- 
dependent of Ai] c , for pseudorapidity separations larger than Arj c « 2 between the 
jet definition cones.Ej This independence of scattering angle in the parton cm. frame 
(for forward scattering) implies that the rapidity gap events follow a 1/t 2 Rutherford 
scattering distribution, just like the non-gap dijet events to which they are normalized 
and which are dominated by single gluon exchange in the t-channel. 

*a 10 3 I 




1 10 100 

^cal 

Fig,_9. Electromagnetic tower multiplicity between widely separated jets as observed by 
DOBil Because of the logarithmic n ca i scale all multiplicities are shifted by +0.5 to the right. 
The excess of low multiplicity (rapidity gap) events above the double negative binomial fit 
(solid and dashed line) is due to color singlet exchange in the t-channel. 

This angular distribution is exactly as would be expected if the gap events were 
produced by t-channel photon exchange. However, the observed rate is too large to be 
consistent with an electroweak origin. The dijet cross section from t-channel photon, Z 
or W exchange is between 10~ 3 and 10~ 4 of the QCD dijet cross section^ in the phase 
space region considered by DO and, hence, cannot account for the one percent fraction 
of dijet events which exhibit a rapidity gap. 

Actually, the discrepancy is even larger because a rapidity gap, even if present at 
the level of a single parton collision, may be covered by double parton scattering (DPS) 
or, more generally, the underlying event. The hardness of the produced jets translates 
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into an impact parameter between the scattering partons which is negligible compared 
to the transverse size of the incident protons. This makes it quite likely that some of 
the other partons in the p and the p will scatter as well, via gluon (color) exchange, 
and the resulting color recombination would lead to the filling of the rapidity gap. Only 
the fraction of dijet events without DPS or an underlying event can be expected to 
preserve a rapidity gap. This fraction, the survival probability P s of rapidity gaps, has 
been estimated by several authorsS™ for both the Tevatron and supercolliders and 
it is expected to be in the 3-30% range. For illustration I will use a value of P s = 10% 
in the following. At the LHC P s can be measured independently, by comparing the 
number of electroweak qq — * qqW or qq — > qqZ events which exhibit a rapidity gap 
with the SM cross sections for these processes.^ 

4-3. Two gluon color singlet exchange and rapidity gaps 

Given that electroweak quark scattering cannot account for the observed rate of 
rapidity gap events at the Tevatron, there must be a QCD source of t-channel color 
singlet exchange. Indeed, the exchange of two gluons in a color singlet state, as depicted 
in Fig. [T(](b), provides a ready source and has been suggested by Low and Nussinov as 
the basic model for the Pomeron.Ell 
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Fig. 10. Feynman graphs for quark scattering via (a) single gluon exchange and (b) exchange 
of two gluons in a color singlet state. 



The rate of rapidity gap events can be understood qualitatively in terms of for- 



ward scattering by single gluon and two gluon color singlet exchange as shown in Fig. |10 
The observed rate was, in fact, predicted by Bjorken.113 For forward qq scattering the 
matrix elements for the dominant helicity configurations are given by 

87ca s s 
Q 2 211 4 * 3 

for single gluon exchange, while color singlet two gluon exchange is dominated by the 
imaginary part of the amplitude in the forward region, 

M singlet « i ^log(Q 2 i? 2 ) ((T b T% 2tl (T b T a ) tm ) smgi£t 

- S7RU ' S as -\ og (Q 2 R 2 ) ti(T b T a )^ll ti(T b T a )^ . (67) 



Q 2 2 ' y ' 3 
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Here R is a cutoff parameter, of the order of the proton radius, which is needed because 
the intermediate phase space integral for the Feynman graph in Fig. 0(b) is infrared 
divergent. The fraction of two gluon color singlet exchange events is then given by 



^singlet 



\\xx\M 



singlet \ 



2\M t 



-log(Q 2 P 2 



hog(Q 2 R 2 ) a s (Q 2 ) 2 



12tt 



(33 - 2n f ) log 



Ji QCD 



(68) 



where we have used the LO formula for the running coupling constant at the scale Q 2 
which is the physical scale of the process. For sufficiently large momentum transfer 
Q 2 = \t\ the two logs are approximately equal and one finds a constant color singlet 
exchange fraction, 



^singlet 



<7n 



An 



33 - 2n 



0.15 



(69) 



for rif = 5. Thus about 15% of all quark-quark scattering events are due to color singlet 
two gluon exchange and may lead to a rapidity gap, provided that no underlying event 
is present. The probability for the latter is given by the survival probability P s and thus 
the expected fraction of dijet events from quark-quark scattering with an observable 
rapidity gap is expected to be§ 



da, 



qq,gap 



/dt 



qq 



da, 



qq, dijet 



/dt 



0.15 P 



0.015 



(70) 



which agrees well with the DO value of 0.0107. This agreement, however, is at least 
partially by accident. In the forward scattering region, the imaginary parts of the color 
singlet exchange amplitudes in quark-gluon and gluon-gluon scattering are the same 
as the quark-quark scattering amplitude, except for a color factor,0 



ImM 



gg, singlet 
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— Im J\A. qg i singlet 



9\ 2 
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(71) 



and, at small t, the same proportionality is found for the single gluon exchange cross 
sections, 



da. 



Ml 



9 da, 



mi 



9\ 2 da, 



<i<j 



, , (72) 

dt 4 dt V47 dt y ' 

Since the singlet exchange cross section is proportional to (Im M. S ingiet) 2 , the gap 
fractions for the three channels are related by 



9 

fgg ~ ~~^qa 



4 J Jqq > 



(73) 



and, thus, the fraction of color singlet exchange events is substantially larger than 
the 15% expected for quark-quark scattering alone. Using the measured rapidity gap 
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fraction of ~ 1% this in turn indicates a survival probability well below 10% for the 
Tevatron. 

The rate estimates above assume that two gluon color singlet exchange, i.e. the 
exchange of a Low-Nussinov pomeron, does indeed lead to the same kind of gluon 
radiation pattern as t-channel photon exchange and, thus, can produce a rapidity gap 
after hadronization. However, the intermediate phase space integral for the Feynman 
graph of Fig. [H^b) is dominated by the phase space region where one of the exchanged 
gluons is very soft. Thus the color of the other, hard gluon is screened at large distances, 
i.e. the Low-Nussinov pomeron is an extended object with colored constituents. As a 
result, analogy to single photon exchange may be misleading: gluon radiation may 
resolve the internal color structure of the pomeron. In order to answer this question 
we need to find out whether real gluon emission in Low-Nussinov pomeron exchange 
follows the pattern typical for i-channel gluon exchange or photon exchange.^ Only 
in the second case can we expect Low-Nussinov pomeron exchange to lead to rapidity 
gaps. 

We already have studied the general color structure of the process q^Qi 3 — > 
qi 2 Qi 4 g a - The amplitude contains two color octet pieces Mi and M 2 and the color 
singlet exchange amplitudes M 12 and M 34 . Let us concentrate on the last one, which 
describes color singlet exchange as seen by quark Q. Even for ^-channel gluon exchange 
does this color singlet amplitude exist. However, the rapidity distribution of the emitted 
gluon is markedly different from photon exchange. 

In the QED case the M34 amplitude corresponds to emission of the final state 
gluon off the quark q (A^ and in Fig. [7]). In forward scattering (j] q = +3 in 
Fig. ||(b)) the gluon is radiated between the initial and final state q directions, i.e. at 
Vg ~ Vg- The color %i of the initial quark q is thus transferred to a low mass color triplet 
object which emerges close to the beam direction. At lowest order this is the final state 
q, at 0(a s ) it is the qg system. The situation is thus stable against gluon emission at 
even higher order for the QED case and gluon radiation is suppressed in the rapidity 
range between the two final state quarks. 




Fig. 11. Two of the 31 Feynman graphs contributing to the imaginary part of the color 
singlet exchange amplitude M34. 



In the QCD case M 34 corresponds to emission of the gluon from the quark Q. 
The gluon is preferentially emitted between the initial and the final Q directions, at 
Vg ~ Vq (dash-dotted line in Fig. ||(b)). Thus the color triplet qg system, into which the 
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initial quark q evolves, consists of a widely separated quark and gluon. Higher order 
corrections will lead to strong gluon radiation into the angular region between the two 
and thus also into the rapidity range between the two final state quarks. 

These typical patterns found for t-channel color singlet and color octet exchange 
can now be used as a gauge for the radiation pattern produced in qQ — > qQg scattering 
via the exchange of two gluons in a color singlet state. In the lowest order process. 
qQ — > qQ, the color singlet exchange amplitude is dominated by its imaginary part.c3 
Hence, we may estimate the radiation pattern by calculating the imaginary part of 



the gluon emission amplitude M 34 only. Typical Feynman graphs are shown in Fig. [I I 
Details of the calculation are given in Ref. 35. 

For massless internal gluon propagators the phase space integrals over the qQ, 
qg, and gQ intermediate states are divergent. They can be regularized by replacing 
the massless gluon propagator by a version which avoids unphysical gluon propagation 
over long distances.^ QCD Pomeron models of this kind have been found to give a 
good description of available data.ca These refinements are approximated by using an 
effective gluon mass of m r = 300 MeV in the calculation. 

A second problem arises because some of the contributions to ImM 34 correspond 
to q — > g splitting and subsequent gQ — > gQ scattering via pomeron exchange. These 
contributions cannot be expected to be suppressed when the gluon is emitted between 
the q and the Q directions and thus would mask the radiation off pomeron exchange in 
qQ scattering. These splitting contributions have been subtracted in Ref. 35 to yield 
the square of the pomeron exchange radiation pattern, |ImM34 m | 2 , which is shown in 



Fig. 12 



For high transverse momentum of the emitted gluon (of order of the quark mo- 



menta, see Fig. |12|(b)) the radiation pattern is quite similar to the one obtained for 
single gluon exchange. Hard emitted gluons have too short a wavelength to see the 
screening of the color charge of the harder exchanged gluon by the second, typically 
very soft, exchanged gluon. The Low-Nussinov pomeron thus reveals itself as an ex- 
tended object. Hard gluon emission is able to resolve the internal color structure of the 
pomeron. As the transverse momentum of the emitted gluon is decreased, a qualitative 
transition occurs, as is apparent by comparing the pr g = 2 GeV and 15 GeV cases in 



Fig. [T3. The gluon radiation has too long a wavelength to resolve the internal color 



structure and hence the pomeron appears as a color singlet object. 

In dijet production at the Tevatron, in the kinematic range studied by DO, gluon 
radiation will be dominated by the soft region, with transverse gluon momenta in the 
few GeV range or below, as shown in Fig. [T2|(a). The emission of these soft gluons follows 
a pattern very similar to the one observed for t-channel photon exchange, which leads 
to the formation of rapidity gaps. One concludes that two gluon color singlet exchange 
in dijet events does indeed lead to the formation of rapidity gap events as observed at 
the Tevatron. 

5. Weak Boson Scattering 

One of the major reasons for building a hadron supercollider is the study of weak 
boson interactions in the TeV range. Within the forseeable future this task will have 
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Fig. 12. Rapidity distribution of emitted gluons in uc — > ucg scattering via two gluon color 
singlet exchange as seen by the charm quark. The phase space parameters for the quarks are 
the same as in Fig. [8| and results are shown for a) the case of a soft gluon (prg = 2 GeV) 
and b) a hard gluon (pTg = 15 GeV). For comparison tree level results are shown for gluon 
(dash-dotted lines) and photon exchange (dashed lines). 



to be performed in pp collisions at 14 TeV at the Large Hadron Collider (LHC) at 
CERN. Because of its limited energy the study of weak boson scattering will be a very 
demanding task at the LHC, requiring a full arsenal of tools to isolate qq — > qqVV 
weak boson scattering events from the backgrounds. For example, jets will have to be 
used as a tagging device (as in the top quark search at the Tevatron) and rapidity gap 
techniques, i.e. the exploitation of the t-channel color singlet exchange in the signal, 
should prove useful as well, thus combining the ideas discussed in the previous Sections. 
Before considering the benefits of these techniques in weak boson scattering and the 
search for a heavy Higgs at the LHC, let us start with a brief note on the physics of a 
strongly interacting Higgs sector. 

5.1. Higgs production and longitudinal weak boson scattering^ 

In a more general context, the search for the Higgs boson is part of the quest 
for the basic interactions which are responsible for the spontaneous breaking of the 
SU(2) x U(l) gauge symmetry. The existence of this local symmetry is evidenced by 
the observation of W and Z bosons and the fact that their experimentally determined 
couplings to quarks and leptons agree with the gauge theory predictions. The appear- 
ance of gauge boson masses then requires the gauge symmetry to be broken sponta- 
neously, with some order parameter $ acquiring a vacuum expectation value (v.e.v.). 
The fact that the W to Z mass ratio leads to a p-parameter very close to unity finally 
tells us that this order parameter must transform essentially as an S'L r (2)-doublet. 

Beyond these basic facts the precise nature of the order parameter and of the 
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Fig. 13. Feynman graphs for the two dominant Higgs production processes at the LHC, (a) 
gluon-gluon fusion via a top-quark loop and (b) weak boson fusion. 



interactions which drive its v.e.v. must be determined experimentally. Within the SM 
the order parameter is simply the scalar Higgs doublet field, 




(74) 



where v is the Higgs v.e.v., x an d X° are the Goldstone bosons which are closely 
related to the longitudinal degrees of freedom of the W and the Z, and H is the field 
describing the physical Higgs boson which we would like to detect experimentally. 

Without spontaneous breakdown of the symmetry, SU(2) x U{1) gauge invariance 
forbids any mass terms for the known quarks and leptons as well as for the gauge 
bosons. Thus these masses must be proportional to the v.e.v. Since v appears only in 
the combination v + H in the SM (see Eq. (|74]) ), this implies couplings of the fermions 
and gauge bosons to the Higgs which are proportional to the masses of the former. 
Let us consider the most important of these, the couplings of the Higgs to the gauge 
bosons and to the top quark. They are derived from the kinetic energy term for the 
Higgs doublet field and from the Yukawa coupling of the Higgs field to the top-bottom 
doublet, T = (t^, b L ) T , and the right handed top quark field, t R , 
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(76) 



The interaction terms of Eq. (|75|) imply that the Higgs boson decays predom- 
inantly into the heaviest particles available. For the heavy Higgs scenario, which is 
considered in the following, these are the decay modes H — > W + W~ , H — > ZZ, and 
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Fig. 14. Production cross sections for the SM Higgs boson at the LHC. From Ref. 46. 



H — > tt. For a Higgs mass of 800 GeV, for example, the corresponding branching ratios 
are 

B(H -»• W + W~) « 0.59 , B(H -> ZZ) « 0.29 , £(# -> ft) « 0.12 , (77) 

and the total Higgs decay width is expected to be about 290 GeV. In order to exploit 
the dominant decay modes, one needs to search for the Higgs signal in W + W~ and 
ZZ production events at the LHC, with subsequent leptonic decays of at least one of 
the weak bosons. 

Similar to the Higgs decay modes, the dominant production processes involve 
heavy particles. At LHC energy, the Higgs boson is mainly produced via gluon fusion, 
gg — > H, which proceeds via a top quark loop and is depicted in Fig. |l3|(a), and via 
weak boson fusion, WW, ZZ — > H, which, more precisely, is the electroweak process 



qq — > qqH (see Fig. [15(b)). The corresponding production cross sections are shown in 
Fig. 0. 

While the gluon fusion cross section ranges from 30 to 0.5 pb for 100 GeV < 
run < 800 GeV, weak boson fusion rates are typically smaller by a factor 2.5 to 10, 
the difference becoming smallest at large Higgs boson masses. This somewhat smaller 
cross section, however, is compensated by a number of features which are special to 
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(2) (3) 

Fig. 15. Feynman graphs for the weak boson scattering process ZZ — ► W + W~ 



the weak boson scattering process. First of all the additional two final state quarks will 
frequently manifest themselves as hadronic jets which allows to suppress backgrounds 
by jet tagging. Second the two production processes probe different interactions of 
the Higgs boson. The gluon fusion cross section is proportional to \ 2 , i.e. it measures 
the Yukawa coupling of the produced scalar to the top quark. This Yukawa coupling, 
however, does not identify the produced particle as being connected to the symmetry 
breaking mechanism, it is possible for a generic scalar particle. The HVV coupling, on 
the other hand, at tree level, is only possible if H is indeed the Higgs boson. The weak 
boson fusion rate directly measures the HVV coupling (as does the H — ► VV decay 
width) and thus one would like to measure both gluon and weak boson fusion cross 
sections in order to disentangle the various couplings of H in Eq. ([75]) . 

When considering the full processes qq — > qqWW or qq — > qqZZ, the resonance 



approximation, i.e. considering graph (b) in Fig. 13 only, is only valid for a Higgs width 



Y H « m H . For a heavy (and wide) Higgs resonance the VV — > H{— ► VV) subprocess 



in Fig. |13~Kb) must be replaced by the full weak boson scattering amplitude, as depicted 



in Fig. 15 for the process ZZ —*_fflW. In addition, weak boson bremsstrahlung off the 
quark lines must be considered.^'" 



Very strong cancellations between the various Feynman graphs of Fig. [15] appear 
when all the external vector bosons are longitudinally polarizedJ-S The isospin and 
angular momentum zero contributions from the first two graphs and from the third 
graph individually grow as s 2 jm^. However, the leading terms cancel when the sum 
of the three graphs is considered. Even the s/m^ growth of their sum, however, would 
lead to a violation of partial wave unitarity at a yl of about 1 TeV. At high energies, 
above s = m 2 H , it is the destructive interference with the fourth graph, s-channel Higgs 
exchange, which leads to an acceptable J = I = partial wave amplitude. Thus the 
existence of the Higgs boson, or some other, additional contribution to the weak boson 



scattering amplitude, beyond the vector boson self-coupling graphs (l)-(3) in Fig. [15 
is required by unitarity of the S'-matrix.ca This additional contribution does not have to 
be the s-channel exchange of a scalar object. It might be the t- or w-channel exchange 
of a vector particle which resides in a weak isospin triplet, like in technicolor models. 
Or nature may have realized the unitarization of the weak boson scattering amplitudes 
in yet another manner.0 

Depending on which path is realized in nature, the various weak boson scattering 
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processes like W + W~ -> W + W~, WZ -> WZ, -> W+W+ etc. will exhibit 

markedly different magnitude and shape (energy dependence) of the cross sections. If 
a technirho existed, the WZ channel, for example, would exhibit a resonance at the 
technirho mass, while this channel has a small cross section if a J = I = scalar is 
responsible for the unitarization of the S'-matrix (as in the SM). Thus the study of 
all weak boson scattering channels is important in order to fully probe the symmetry 
breaking sector. In the following I will use the search for a SM heavy Higgs to discuss 
techniques of background suppression. One should keep in mind, however, that these 
methods are more general and can be applied to study any of the weak boson scattering 
processes. 

5.2. Forward jet-tagging 

A characteristic feature of weak boson scattering events are the two accompanying 
quarks (or antiquarks) from which the "incoming" Ws or Zs have been radiated (see 
Fig. [T3|(b)). In general these scattered quarks will give rise to hadronic jets. By tagging 
them, i.e. by requiring that they are observed in the detector, one hopes to obtain a 
powerful background rejection tool£3'c3 Whether such an approach can be successful 
depends on the properties of the tagging jets: their typical transverse momenta, their 
energies, and their angular distributions. 

Similar to the emission of virtual photons from a high energy electron beam, the 
incoming weak bosons tend to carry a small fraction of the incoming parton energy .ca 
At the same time the incoming weak bosons must carry substantial energy, of order 
m#/2 = myy/2, in order to produce a weak boson pair of large invariant mass. Thus 
the final state quarks in qq — > qqVV events will carry very high energies, of order 
1 TeV or even higher. This is to be contrasted with their transverse momenta, which 
are of order p T « m w . This low scale arises because the weak boson propagators in 
Fig. 0(b) introduce a factor 

Dv{q 2 ) = 2 ~ l 2 « 2 ] 2 (78) 
q z — my pj, + m v 

into the production amplitudes and suppress the qq — > qqH cross section for quark 
transverse momenta above my. The modest transverse momentum and high energy of 
the scattered quark corresponds to a small scattering angle, typically in the 1.5 < r\ < 
4.5 pseudorapidity region. 

These general arguments are confirmed by Fig. [16], where the transverse momen- 
tum and pseudorapidity distributions of the two potential tagging jets are shown for 
the production of a m H = 800 GeV Higgs boson at the LHC. One finds that one 
of the two quark jets has substantially lower median (« 30 GeV) than the other 
(~ 80 GeV). As a result, double jet-tagging-i.e^the requirement that both scattered 
quarks are visible as jets, proves quite costly£jtL3 unless jets with transverse momenta 
around 30 GeV can be identified in the forward region. A different approach is single 
forward jet-tagging which relies only on the higher pt tagging^jet and thus proves more 
effective for higher transverse momentum thresholds." 
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Fig. 16. Transverse momentum and pseudorapidity distributions of the two (anti)quark jets 
in qq — > qqH events at LHC energy. Shown are (a) da/dprj for the highest (dashed curve) 
and lowest pt jet (solid curve) and (b) da/d\r]j\ for the most forward (dashed curve) and the 
most central jet (solid curve). 



This transverse momentum threshold needs to be set such that the probability for 
seeing a fake tagging-jet becomes small. At full LHC luminosity of £ = 10 34 cm~ 2 sec _1 
one of the most important sources of such background jets are additional pp — > jjX 
events in the same bunch crossing as the pp — > VVX signal event. With a bunch cross- 
ing happening every 25 nsec, a process with a cross section of a e ff = [C 25nsec] _1 = 
4 mb will, on average, happen during each bunch crossing. Using Poisson statistics, 
the probability to get a random jet of transverse momentum p?j > Pr,cut m any event 
previously selected is then given by 

P,W) = 1 - exp L*»fa « -n(PTi>PT,^) (7g) 

V Veff J CTeff 

With single-jet cross sections of about 0.5-1 mb above pTj = 20 GeV in the relevant 
rapidity range, overlapping events in a single bunch crossing may produce fake tagging- 
jets above 20 GeV px with about 15-25% probability. This estimate agrees with the 
results of a more detailed analysis of overlapping events at the LHC.il 

In the following let us sidestep this question of how low a transverse momentum 
threshold can be achieved in forward jet-tagging. By restricting ourselves to single 
forward jet-tagging the jet px threshold can be set high enough to avoid the fake 
jet problem from overlapping events. As an example we will consider the search for 
H — > W + W~ — > i + vl~v decays as discussed in Ref. 18. This particular decay mode 
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suffers from severe backgrounds due to tt — > W + W~bb decays as well as qq — > W + W~ 
events (referred to as WW QCD background in the following) and thus is well suited 
to assess the power of jet-tagging and minijet veto techniques. 

Since we are interested in the decay of a very heavy Higgs boson, the two charged 
W decay leptons will emerge with high transverse momentum, in the central region of 
the detector, and they will be well isolated from additional jets. Thus we require the 
presence of two charged leptons (£ = e, fj) with 

p T£ >50GeV, M<2, R tj = ^ - Vj ) 2 + (0, - 0,) 2 > 0.7 . (80) 

Here pxt denotes the lepton transverse momentum and rjg is its pseudorapidity. The 
Rfj > 0.7 separation cut forbids a parton (jet) of p? > 20 GeV in a cone of radius 0.7 
around the lepton direction. 

The lepton pt cut in Eq. QgOQ is not in itself sufficient to focus on the production 
of two Ws of large transverse momenta and large W^-pair invariant mass. A variable 
which helps to substantially suppress W bremsstrahlung backgrounds is Apree, the 
difference of the charged lepton transverse momentum vectors£il We thus require 

Ap Tee = \p Th - p Ti2 1 > 300 GeV , m u > 200 GeV . (81) 

The additional cut on the dilepton invariant mass removes possible backgrounds from 
Z leptonic decays. It is largely superceded by the the Apxee cut, however. 

Cross sections for events satisfying the lepton acceptance criteria of Eqs. ( PU| , PT]) 
are listed in the first column of Table |2| for the case of a ra# = 800 GeV Higgs 
boson and the qq — > W + W~ and tt production backgrounds. In the calculation of the 
qq — > qqWW signal all W bremsstrahlung graphs have been includedEZI'El since the 
narrow Higgs width approximation is not appropriate any more for the large masses 
considered here. These additional graphs introduce an "electroweak background" which 
is the same as in the case of a light Higgs boson. Since a light Higgs cannot decay via 
H — > WW, we can use the SM qq — > qqWW cross section for, say, mu = 100 GeV as 
an estimate of this electroweak background. The heavy Higgs signal is then defined as 
Basic = Ba(m H ) - Ba(m H = 100 GeV). 

Without any jet-tagging the backgrounds are overwhelming: the top quark back- 
ground alone exceeds the signal by a factor of 300. How much does this situation 
improve by forward jet-tagging? A large fraction of the signal (55%) is retained by 
requiring the existence of a very energetic, forward jet of moderate p?, in the phase 
space region 

Ef g > 500 GeV , 1.5 < |rf g | < 4.5 , p^f > 50 GeV . (82) 

In addition the tagging-jet should be well separated from the W decay leptons, 

min |?7* as — %| > 1.7 , (83) 

Since we are now requiring additional jet activity, the previous background cal- 
culations are not sufficient any more. Rather 0(a s ) real emission corrections must be 
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Table 2. Signal and background cross sections Ba (in fb) for the Higgs search in qq — ► 
qqWW, WW — > £z^i/ events. Results are shown for increasingly stringent cuts. From Ref. 18. 



lepton cuts only + tagging jet 



WW(jj) 
ti(jj) 

m H = 100 GeV 
m H = 800 GeV 



[Eq. m-(M)} [Eq- (H)] 



27.4 
640 
1.18 
3.4 



1.73 
57 
0.56 
1.79 



+ lepton- 
tagging jet 
separation 

[Eq. (M 
0.57 
25 
0.29 
1.31 



+ minijet veto 
(pr.veto = 
20 GeV) 

[Eq- (jj)] 
0.13 
0.47 
0.18 
0.97 



signal: 

m H = 600 GeV 
m H = 800 GeV 
mu = 1 TeV 



2.2 



1.23 



1.02 



0.78 
0.79 
0.62 



considered, i.e. ttj production and processes like qq — > W + W~g and qg — > W + W~q 
need to be considered at this level. The signal and background cross sections after 
the cuts of Eqs. fl8"2p and ( |S3"D are listed in the second and third columns of Table ^], 
respectively. Single forward jet-tagging has reduced the backgrounds by a factor 25- 
50, while keeping 45% of the signal. This reduction factor is almost sufficient for the 
WW QCD background, but an additional large suppression factor is needed for the 
top quark decays. Fortunately, the 6-quarks in t — > bW decay frequently give rise to 
additional jets and a veto on them,@ above pxb = 25 GeV, will substantially reduce this 
background.!^ One needs to be careful, however, since at such low transverse momenta 
the production of minijets via the emission of additional gluons cannot be neglected at 
the LHC. 

5.3. Rapidity gaps and minijet veto 

The most rigorous veto on 6-quarks from top decays would be a rapidity gap 
requirementJll0 as discussed for dijet events at the Tevatron in Section [4.2| . A rapidity 
gap trigger, i.e. the selection of events without any hadronic activity in some rapid- 
ity range around the two charged leptons which arise from the two Ws, would make 
maximal use of the different color structure of signal and background processes. In a 
weak boson scattering event no color is exchanged between the initial state quarks. 
Color coherence between initial and final state gluon bremsstrahlung then leads to 
a suppression of hadron production in the central region, between the two scattered 
quarks.EJUZlEa The tt production and WW QCD backgrounds, on the other hand, in- 
volve color exchange between the incident partons and, as a result, gluon radiation into 
the central region dominates. This gluon radiation then results in substantial hadronic 
activity in the vicinity of the W decay leptons. 
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While a rapidity gap trigger should suppress the backgrounds well below the 
signal level, it is not practical at the LHC, due to the small signal rate and small 
survival probability of the signal. In order for a rapidity gap to be visible no additional 
scattering process may occur in the detector at the same time, neither in the same pp 
collision which produces the Higgs boson, nor by other proton pairs which collide in the 
same bunch crossing. The first condition is parametrized by the survival probability, 
P s , which was discussed in Section fO| , and which we expect to be below 10% at 
the LHC. 081 The second condition implies that one must run the LHC well below 
design luminosity, at an average of one pp scattering event per bunch crossing, or about 
C = 5 • 10 32 cm _2 sec _1 . Even at this low luminosity only P\ = e -1 = 37% of all bunch 
crossings have exactly one pp scattering event, leading to an effective yearly integrated 
luminosity of about 2 fb _1 . Given the signal cross section of 1.02 fb in the third column 
of Table |2] and an (optimistic) survival probability of P s = 10%, one would register 
only one Higgs signal event every five years! Note that the cuts imposed on the leptons 
and tagging-jet arenot responsible for this negative result: 22% of all signal events are 
accepted by them.Ha It is the combination of small survival probability and additional 
effective luminosity reduction by a factor 50 which kills the signal. Given the small 
weak boson scattering cross section, any rapidity gap strategy at the LHC must raise 
this effective survival probability from 1/500 to something of order unity, i.e. it must 
work at full design luminosity and must tolerate the presence of an underlying event. 

The way out may be to use minijets rather than soft hadrons to define the rapid- 
ity gap. Beyond the different angular distribution of gluons in signal and background 
events, a second distinction is the momentum scale Q of the hard process which gov- 
erns the radiation. In longitudinal weak boson scattering the color charges, carried by 
the incident quarks, receive momentum transfers given by the transverse momenta of 
the final state quarks, which typically are in the Q = 30 to 80 GeV range. For the 
background processes, on the other hand, the color charges receive a much larger mo- 
mentum kick, of the order of the weak boson pair mass or even the parton center of 
mass energy of the event, i.e. Q ~ 1 TeV. Extra parton emission is suppressed by a 
factor f s = a s \n (Q 2 jp\ min ) , where pr,min is the minimal transverse momentum re- 
quired for a parton to qualify as a jet. The jet transverse momentum scale below which 
multiple minijet emission must be expected is set by f s = 0(1), and this scale is in the 
30-50 GeV range for the backgrounds but well below 10 GeV for the signal. 

For a more precise determination of the typical minijet pt scale we must carry 
signal and background calculations still one order (in a s ) higher, i.e. the emission of 
one additional parton must be included. Thus we are led to consider qq — > qqWWg 
and crossing related _processesEj for the weak boson scattering signal, and for the back- 



grounds ttjj events^ and WWjj events™ must be simulated. Again the methods 
discussed in Section || are very useful in performing these tasks. The results of these 
calculations confirm the qualitative arguments made above. For the hard lepton and 
tagging-jet cuts of Eqs. (|80|)-(|83D, the WWjj cross section, for example, is equal to the 
lowest order WWj cross section for a minimal parton transverse momentum of 37 GeV. 
Thus multiple minijets in the 20-50 GeV range are expected for the WW QCD back- 
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ground. Similar results are found for the tt background while this "saturation scale" 
for the signal is well below 10 GeV. 

For a quantitative estimate of the minijet emission probability, we must use the 
higher order programs (which include emission of soft partons) in regions of phase 
space where the n + 1 jet cross section saturates the rate for the hard process with n 
jets. As the pr of the softest jet is lowered to values where a{n + 1 jet) ~ cr(n jet), 
fixed order perturbation theory breaks down and multiple soft gluon emission (with 
resummation of collinear singularities into quark and gluon structure functions, etc.) 
needs to be considered in a full treatment. For the complex processes considered here 
such calculations are not yet possible. Instead one can employ the so called "truncated 
shower approximation" (TSA) to normalize the higher order emission calculations.^ 
The tree-level n + 1 jet differential cross section da(n + 1 j)tl is replaced by 

da(n + 1 j ) tsa = da(n + 1 j) TL (l - e'^-^^s^ , (84) 

with the parameter ptsa properly chosen to correctly reproduce the lower order n jet 
cross section when integrated over a given phase space region of this hard process. Here 
PTj min is the smallest transverse momentum of any of the final state massless partons. 
As PTj,min — ► the final factor in Eq. (|84| ) acts as a regulator, which allows to integrate 
the n + 1 parton cross section over the full phase space region and simulate n-jet and 
n + 1-jet events simultaneously. 

The results of such a calculation^ are depicted in Fig. |T7[ The angular distribution 
of the jet (parton with p? > 20 GeV) which is closest to the leptons (more precisely, 
closest to the average lepton rapidity fj = (r] e + + %-)/2) is shown in Fig. [13(a). The 
background processes favor emission close to the leptons. For the Higgs signal this 
closest jet is typically the second quark in the qq — > qqH process and not soft gluon 
radiation. It is this different angular distribution which is at the heart of the rapidity 
gap trigger. 

A good compromise between strong background rejection and high signal accep- 
tance is achieved by vetoing jets in the veto region defined by 

P v tT > W>eto , rj] eto e [rjf n - 1.7, rg*] or faf g , r,?^ + 1.7] . (85) 

The veto probability as a function of the cut value pr,veto is shown in Fig. [H](b). Even 
though these results were obtained in the truncated shower approximation and a more 
precise modeling is needed, they clearly demonstrate that, in the central region, the 
backgrounds have a much higher probability to produce additional minijets from QCD 
radiation than the weak boson scattering signal. Even for the tt background, where a 
strong suppression is obtained by vetoing the central 6-quark jets arising from the top 
decays, the veto on the jet activity from soft QCD radiation provides an additional 
suppression by a factor 2 for pT, ve to = 20 GeV. Final cross section values for signal and 
backgrounds are given in the last column of Table §. Should minijet vetoing be possible 
at the LHC for even smaller pT, ve to values then the tt background to H — > WW events 
can effectively be eliminated. 
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Fig. 17. Rapidity and transverse momentum distributions of secondary jets. In (a) Ad- 
measures the pseudorapidity distance of the jet closest to the leptons from the average lepton 
rapidity fj. Negative values of Ar/ij correspond to soft jets on the opposite side of the leptons 
with respect to the tagging jet. The dashed line shows the distribution for the electroweak 
background as defined by the mn = 100 GeV case. The ttjj background has been scaled down 
by a factor 30. The probability to find a veto jet candidate above a transverse momentum 
PT,veto) i n the veto region of Eq. (|8q) , is shown in (b). From Ref. 18. 



6. Conclusions 

Rapidity gaps in hard dijet events at the Tevatron and minijet patterns in weak 
boson scattering at the LHC have their common origin in processes which are dom- 
inated by t-channel color singlet exchange. Both can be understood in terms of soft 
gluon radiation patterns which are determined by the color structure of the underlying 
hard scattering event. At the lower energy of the Tevatron the soft gluon radiation 
pattern is reflected in the distribution of soft hadrons in events without a residual min- 
imum bias event (underlying event). Going up one order of magnitude to LHC energy, 
"soft gluons" can be hard enough to be seen as distinct minijets, and this opens new 
strategies for the study of weak boson scattering events. 

The theoretical description of these multi-parton processes has become possible 
with the advent of amplitude techniques. By directly evaluating polarization ampli- 
tudes numerically, fast numerical programs can be designed which remain efficient for 
processes with large numbers of Feynman graphs. In addition, code generation for SM 
processes has now been automated and this allows us to concentrate on the physics 
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issues in the analysis of multi-parton processes. 

Via the study of rapidity gap events we are obtaining important information for 
future LHC physics studies already now at the Tevatron. Similarly, multiple minijet 
emission, which is expected to be quite common at the LHC, can already now be 
studied in the hardest "dijet" events at the Tevatron.HI! Both examples show the 
important interplay between present experiments and new theoretical ideas and tools 
in multi-parton analysis. Some of these developments were motivated by the search 
for new physics at future colliders. At the same time they aid in the preparation of 
experiments at these machines. 

Acknowledgements 

Special thanks go to the organizers for making this most enjoyable and stimulating 
meeting possible. In preparing this talk I have benefitted from many discussions with 
V. Barger, J. D. Bjorken, A. Brandt, F. Halzen, W. Long, R. Phillips, and D. Summers 
and I want to thank them for sharing their insights with me. I would like to thank 
T. Stelzer for teaching me about MadGraph and for a critical reading of Section 3. 
This research was supported in part by the University of Wisconsin Research Com- 
mittee with funds granted by the Wisconsin Alumni Research Foundation, and by the 
U. S. Department of Energy under Grant No. DE-FG02-95ER40896. 



37 



References 

1. UA1 Collaboration, G. Arnison et al, Phys. Lett. 122B (1983) 103; UA2 Collabo- 
ration, M. Banner et al, Phys. Lett. 122B (1983) 476. 

2. CDF Collaboration, F. Abe et al., Phys. Rev. D50 (1994) 2966; Phys. Rev. D52 
(1995) 2605. 

3. CDF Collaboration, F. Abe et al., Phys. Rev. Lett. 74 (1995) 2626. 

4. DO Collaboration, S. Abachi et al, Phys. Rev. Lett. 74 (1995) 2422; Phys. Rev. 
Lett. 74 (1995) 2632; Phys. Rev. D52 (1995) 4877. 

5. For a recent review see e.g. X. Tata, lectures given at TASI 1995, Univ. of Hawaii 
preprint UH-511-833-95 (1995), |hep-ph/ 95 10287 . 

6. P. De Causmaeker, R. Gastmans, W. Troost, and T. T. Wu, Nucl. Phys. B206 
(1982) 53; F. A. Berends, R. Kleiss, P. De Causmaeker, R. Gastmans, W. Troost, 
and T. T. Wu, Nucl. Phys. B206 (1982) 61; B264 (1986) 265 and references therein; 

D. Danckaert, P. De Causmaeker, R. Gastmans, W. Troost, and T. T. Wu, Phys. 
Lett. 114B (1982) 203. 

7. Z. Xu, Da-Hua Zhang, and L. Chang, Tsinghua University Preprint, Beijing, TUTP- 
84/4, TUTP-84/5, TUTP-84/6, and Nucl. Phys. B291 (1987) 392; J. F. Gunion 
and Z. Kunszt, Phys. Lett. 161B (1985) 333; R. Kleiss and W. J. Stirling, Nucl. 
Phys. B262 (1985) 167. 

8. K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B274 (1986) 1; ibid. B313 (1989) 
560. 

9. H. Murayama, I. Watanabe, and K. Hagiwara, KEK Report 91-11 (1992). 

10. T. Stelzer and W. F. Long, Comput. Phys. Commun. 81 (1994) 357. 

11. For a recent review see e.g. J. Bagger et al., Phys. Rev. D49 (1994) 1246; Phys. 
Rev. D52 (1995) 3878. 

12. R. N. Cahn et al., Phys. Rev. D35 (1987) 1626; V. Barger et al., Phys. Rev. D37 
(1988) 2005; R. Kleiss and W. J. Stirling, Phys. Lett. 200B (1988) 193. 

13. U. Baur and E. W. N. Glover, Nucl. Phys. B347 (1990) 12; U. Baur and 

E. W. N. Glover, Phys. Lett. B252 (1990) 683; D. Froideveaux, in Proceedings 
of the ECFA Large Hadron Collider Workshop, Aachen, Germany, 1990, edited by 
G. Jarlskog and D. Rein (CERN report 90-10, Geneva, Switzerland, 1990), Vol II, 
p. 444; M. H. Seymour, ibid, p. 557; V. Barger, K. Cheung, T. Han, J. Ohnemus, 
and D. Zeppenfeld, Phys. Rev. D44 (1991) 1426. 

14. V. Barger, K. Cheung, T. Han, and R. J. N. Phillips, Phys. Rev. D42 (1990) 3052. 

15. V. Barger, K. Cheung, T. Han, and D. Zeppenfeld, Phys. Rev. D44 (1991) 2701; 
D48 (1993) 5444E; D48 (1993) 5433. 

16. Y. L. Dokshitzer, V. A. Khoze, and S. Troyan, in Proceedings of the 6th Interna- 
tional Conference on Physics in Collisions, (1986) ed. M. Derrick (World Scientific, 
Singapore, 1987) p. 365. 

17. J. D. Bjorken, Int. J. Mod. Phys. A7 (1992) 4189; Phys. Rev. D47 (1993) 101; 
preprint SLAC-PUB-5823 (1992). 

18. V. Barger, R. J. N. Phillips, and D. Zeppenfeld, Phys. Lett. B346 (1995) 106. 

19. DO Collaboration, S. Abachi et al., Phys. Rev. Lett. 72 (1994) 2332. 

20. CDF Collaboration, F. Abe et al., Phys. Rev. Lett. 74 (1995) 855. 



38 



21. DO Collaboration, S. Abachi et al., Phys. Rev. Lett. 76 (1996) 734. 

22. P. Nason, S. Dawson, and R. K. Ellis, Nucl.Phys. B303 (1988) 607; E. Laenen, 
J. Smith, and W. L. van Neerven, Nucl.Phys. B369 (1992) 543; Phys. Lett. B321 
(1994) 254. 

23. S. Catani, M. Mangano, P. Nason, and L.Trentadue, preprint CERN-TH-96-21 
(1996), |Ee"p-ph/9602208| . 

24. V. Barger, T. Han, J. Ohnemus, and D. Zeppenfeld, Phys. Rev. D40 (1989) 2888; 
F. A. Berends, W. T. Giele, and H. Kuijf, Nucl. Phys. B321 (1989) 39; F. A. Berends 
et al., Phys. Lett. B224 (1989) 237. 

25. F. A. Berends, H. Kuijf, B. Tausk, and W. T. Giele, Nucl. Phys. B357 (1991) 32. 

26. K. Kleinknecht, Talk presented at the XXXth Rencontre de Moriond "Electroweak 
interactions and unified theories", Les Arcs, Savoie, 11-18 March 1995. 

27. G. Peter Lepage, J. Comput. Phys. 27, 192 (1978). 

28. S. Kawabata, Comput. Phys. Commun. 41 (1986) 127. 

29. U. Baur and D. Zeppenfeld, Phys. Rev. Lett. 75 (1995) 1002; E. N. Argyres et al., 
Phys. Lett. B358 (1995) 339. 

30. K. Hagiwara et al., Nucl. Phys. B282 (1987) 253. 

31. J. Kublbeck, M. Bohm, and A. Denner, Comput. Phys. Commun. 60 (1990) 165. 

32. R. Mertig, M. Bohm, and A. Denner, Comput. Phys. Commun. 64 (1991) 345; 
A. Hsieh and E. Yehudai, Comput. Phys. 6 (1992) 253; T. Ishikawa et al, KEK 
Report 92-19 (1993); E. E. Boos et al., Int. J. Mod. Phys. C5 (1994) 615. 

33. V. Barger et al, Phys. Lett. B338 (1994) 336. 

34. D. Zeppenfeld, Int. J. Mod. Phys. A3 (1988) 2175; K. Hagiwara and D. Zeppenfeld, 
Nucl. Phys. B313 (1989) 560; V. Barger, A. L. Stange, and R. J. N. Phillips, Phys. 
Rev. D45 (1992) 1751. 

35. H. Chehime and D. Zeppenfeld, preprint MAD/PH/814 (1994), |hq>ph/9401244 . 

36. Y. L. Dokshitzer et al, Rev. Mod. Phys. 60 (1988) 373, and references therein. 

37. R. S. Fletcher and T. Stelzer, Phys. Rev. D48 (1993) 5162. 

38. H. Chehime et al, Phys. Lett. B286 (1992) 397. 

39. E. Gotsman, E. M. Levin, and U. Maor, Phys. Lett. B309 (1993) 199. 

40. H. Chehime and D. Zeppenfeld, Phys. Rev. D47 (1993) 3898. 

41. F. E. Low, Phys. Rev. D12 (1975) 163; S. Nussinov, Phys. Rev. Lett. 34 (1975) 
1286. 

42. See e.g. J. R. Cudell and B. U. Nguyen, Nucl. Phys. B420 (1994) 669 and references 
therein. 

43. P. V. Landshoff and O. Nachtmann, Z. Phys. C35 (1987) 405; 

44. F. Halzen, G. I. Krein, and A. A. Natale, Phys. Rev. D47 (1992) 295; M. B. Gay 
Ducati, F. Halzen, and A. A. Natale, Phys. Rev. D48 (1993) 2324. 

45. For a more complete discussion see e.g. J. F. Gunion, H. E. Haber, G. Kane, and 
S. Dawson, The Higgs Hunter's Guide, Addison- Wesley (1990). 

46. J. F. Gunion, A. Stange, and S. S. D. Willenbrock, preprint UCD-95-28 (1995), 
hep-ph/ 9602238| . 

47. D. A. Dicus and R. Vega, Phys. Rev. Lett. 57 (1986) 1110; Phys. Rev. D37 (1988) 
2427; J. F. Gunion, J. Kalinowski, and A. Tofighi-Niaki, Phys. Rev. Lett. 57 (1986) 
2351; U. Baur and E. W. N. Glover, Phys. Lett. B252 (1990) 683. 



39 



48. D. A. Dicus and V. S. Mathur, Phys. Rev. D7 (1973) 3111; B. W. Lee, C. Quigg, 
and H. Thacker, Phys. Rev. D16 (1977) 1519; M. Veltman, Acta Phys. Pol. B8 
(1977) 475; M. S. Chanowitz and M. K. Gaillard, Nucl. Phys. B261 (1985) 379. 

49. See e.g. V. Barger and R. Phillips, Collider Physics, Addison- Wesley (1987). 

50. V. Barger et al., Phys. Rev. D46, 2028 (1992). 

51. D. Dicus, J. F. Gunion, and R. Vega, Phys. Lett. B258, 475 (1991); D. Dicus, 
J. F. Gunion, L. H. Orr, and R. Vega, Nucl. Phys. B377, 31 (1991). 

52. G. Ciapetta and A. DiCiaccio, in Proceedings of the ECFA Large Hadron Collider 
Workshop, Aachen, Germany, 1990, edited by G. Jarlskog and D. Rein (CERN 
report 90-10, Geneva, Switzerland, 1990), Vol II, p. 155. 

53. J. F. Gunion et al, Phys. Rev. D40 (1989) 2223. 

54. A. Duff and D. Zeppenfeld, Phys. Rev. D50 (1994) 3204. 

55. A. Stange, private communication. 

56. V. Barger, T. Han, J. Ohnemus, and D. Zeppenfeld, Phys. Rev. D41 (1989) 2782. 

57. V. Barger and R. J. N. Phillips, Phys. Rev. Lett. 55 (1985) 2752; H. Baer, 
V. Barger, H. Goldberg, and R. J. N. Phillips, Phys. Rev. D37 (1988) 3152. 

58. CDF Collaboration, F. Abe et al., Phys. Rev. Lett. 75 (1995) 608. 

59. D. Summers and D. Zeppenfeld, preprint MADPH-95-904 (1995), |hep-ph/ 9509206 . 



40 



